Z 2 invariant protected bound states in topological insulators 
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We present an exact solution of a modifed Dirac equation for topological insulator in the presence 
of a hole or vacancy to demonstrate that vacancies may induce bound states in the band gap of 
topological insulators. They arise due to the Z2 classification of time-reversal invariant insulators, 
thus are also topologically-protected like the edge states in the quantum spin Hall effect and the 
surface states in three-dimensional topological insulators. Coexistence of the in-gap bound states and 
the edge or surface states in topological insulators suggests that imperfections may affect transport 
properties of topological insulators via additional bound states near the system boundary. 

PACS numbers: 73.20.-r, 73.20.Hb, 74.43.-f 



Topological insulators are narrow-band semiconduc- 
tors with band inversion generated by strong spin-orbit 
coupling [lj. They are distinguished from the ordinary 
band insulators according to the Z2 invariant classifica- 
tion of the gapped band insulators due to the time re- 
versal symmetry. The variation of the I2 invariant at 
their boundaries will lead to the topologically protected 
edge or surface states with the gapless Dirac energy 
spectrum [2 Imperfections, such as impurity, vacancy, 
and disorder, are inevitably present in topological insu- 
lators. Owing to the time-reversal symmetry, an exciting 
feature of topological insulator is that its boundary states 
are expected to be topologically protected against weak 
non- magnetic impurities or disorders [8|,l9(. This provoked 
much interest on the single impurity problem on the sur- 
face of a topological insulator, starting with gapless Dirac 
modelfiol- 14 1. However, reminding that the boundary 
state is only a manifestation of the topological nature of 
bulk bands, it should also start with the examination of 
the host bulk to know how the imperfections affect the 
electronic structure. It is well known that single impurity 
or defect can induce bound states in many systems, such 
as in the Yu-Shiba state in s-wave superconductor ll| 16 1 
and in d-wave superconductors [lTj- Topological defects 
were discussed in the B-phase of 3 He superfiuid[l8| and 
topological insulators and superconductors [l9[ . Here we 
report that bound states can form around a single va- 
cancy in the bulk energy gap of topological insulators. 
These bound states are found to have the same origin 
as boundary states due to the Z2 classification, thus are 
also topologically protected. 

The formation of the topological bound states can be 
readily illustrated by reviewing the quantum spin Hall 
effect in two-dimensional (2D) topological insulators [20I— 
22j, in which strong spin-orbit coupling twists the bulk 
conduction and valence bands, leading to a nontrivial Z2 
index. As the Z2 varies across the edge, edge states arise 
in the gap with the gapless Dirac dispersion. Unlike the 
quantum Hall effect in a magnetic field, spin-orbit cou- 
pling preserves the time reversal symmetry, so the result- 
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FIG. 1: Schematic description of the formation of vacancy- 
induced in-gap bound states in two-dimensional topological 
insulators, [(a) and (b)] A pair of helical edge states traveling 
along the edge of a 2D topological insulator with the gapless 
Dirac dispersion, [(c) and (d)] When the edge is bent into 
a hole, the helical edge states evolve to circulate around the 
hole, [(e) and (f)] The circulating edge states may develop 
into topologically-protected bound states as the hole shrinks 
into a point or being replaced by a vacancy. The same physics 
is expected to happen in one and three dimensions. 



ing edge states appear in helical pairs, in which one state 
is the time-reversal counterpart of the other, propagating 
along opposite directions and with opposite spins (Fig. 
QJ. Now imagine that the system edge is rolled into a 
hole, the edge states will circulate around the hole as 
the periodic boundary condition along the propagating 
direction remains unchanged. While shrinking its radius 
of the hole, most of the edge states will be expelled into 
the bulk bands as the energy separation of the states be- 
comes larger and larger, and it is found that at least two 
degenerate pairs of the states will be trapped to form the 
bound states in the gap while the hole will evolve into 
a point defect. This mechanism of the formation of the 
bound states can be realized in topological insulator in 
all the dimensions. 
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We will employ a modified Dirac model to provide a 
unified description of topological insulators in various di- 
mensions 



H = vp ■ a 



\mv 



Bp 2 ) /3. 



(1) 



The modification comes from the quadratic correction in 
momentum —Bp 2 to the band gap mv 2 . pi = —ihdi is the 
momentum operator (i g {x, y, z}), p 2 = p 2 + p 2 + p 2 ,, v 
and m have the dimension of the speed and mass, respec- 
tively. B has the dimension of m . The Dirac matri- 
ces satisfy the anticommutation relations onotj — —a^cti 
(i 7^ j), di/3 = —fiat and a 2 — (3 2 — 1. One represen- 
tation of the Dirac matrices in three spatial dimensions 
can be expressed as a set of 4 x 4 matrices 



a. 



1 0o, 



(2) 



where <7i =x ,y,z are the Pauli matrices, <tq is the 2x2 unit 
matrix, and ® represents the Kronecker product. This 
modified Dirac Hamiltonian preserves the time reversal 
symmetry QHqQ~ 1 = Hq under the time reversal oper- 
ation O = —ia x a z K , where K is the complex conjugate 
operator. This equation has the identical mathematical 
structure as the effective models for the quantum spin 
Hall effect and 3D topological insulator [21, 23-25[. It be- 
comes topologically non-trivial if mB > 0, while topolog- 
ically trivial if mB < according to the Z2 classification 
of insulato rs I2H I27I . Due to the bulk-boundary corre- 
spondence [28l l29j. there always exist topologically pro- 
tected boundary states at the open boundaries of a topo- 
logical insulator, where the Z2 invariant changes from 
nontrivial to trivial. This feature can be well described 
by the modified Dirac model when mB > 0. Starting 
from this modified Dirac model, we are now ready to 
explore existence of the in-gap bound states induced by 
a single vacancy by presenting an exact solution of the 
modified Dirac model. 

In two dimension (p z = 0) the equation can be reduced 
into two independent 2x2 hamiltonians 



Ha 



(mv 2 



Bp 2 )a z 



hv(p x a x ±p y (Ty 



(3) 



with h_ the time-reversal counterpart of h+, which was 
already studied in the Hg Te quantum well[21|, and in 
the thin films of Bi-^Se^ |24l. [25| . It is convenient to adopt 
polar coordinates (x, y) = r(cos<p, sirup) in two dimen- 
sions. Here these equations are solved under the va- 
cancy boundary conditions [Fig. Ufa)], i.e., the center 
of the 2D topological insulator is punched with a hole 
of radius R, thus the wavefunction is required to vanish 
at r = R and r = +00. Due to the rotational sym- 
metry, the z-component of the total angular momentum 
j z = —ihdg + (h/2)<j z provides a good quantum number, 
labeled by a half-integer mj <G {±1/2, ±3/2, ...}, which 
can be used to characterize the bound states. In this way 
the equation is reduced to a set of ID radial equations, 
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FIG. 2: Two-dimensional topologically protected bound 
states, (a) A 2D topological insulator with a hole of radius 
R at the center, (b)-(c), Energies (E in units of the band 
gap A) of in-gap bound states circulating around the hole 
as functions of the hole radius, rrij is the quantum num- 
ber for z-component of total angular momentum of the cir- 
culating bound states. In (b), m — v = B = h = 1 and 
A = \/3 in (c), mv 2 = -10 meV, Bh 2 = -686 meV-nm 2 , and 
hv = 364.5 meV nm, A = 20meV, adopted from Ref. HJ; 
in (d), mv 2 = 0.126 eV, Bh 2 = 21.8 eVA 2 , hv = 0.294 eVA, 
A = 0.252eV, adopted from Ref. H}. 



which can be solved exactly. The trial wave function has 
the form (1^1,1^2) eXr ■ The secular equations of the in- 
determinate coefficients (V'i,V ; 2) T gi ye f° ur roots of A„ 
(— ±Ai, ±A2) in terms of the energy E, 



\ 2 _ 
A l,2 — 



2B 2 h 2 



[1 - 2mB± y/i - AmB + AB 2 E 2 /v 4 }. (4) 



Using the boundary conditions at r — R and r = +00 
we finally arrive at the transcendental equation for the 
bound state energies 

Xl + K mi+ ^ R ) A 2 + ^™ 3+ i(A 2 i?) 
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and the wavefunction ty m . (r, 9) for h+ turns out to have 
the form 
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where K n (x) is the modified Bessel function of second 
kind. 

In Fig. [2] (b)-(d), we show the bound-state energies 
as functions of R for an ideal case [(b), mB = 1], for 
the HgTe quantum well (mB = 0.05) (SU, and for a 2 
quintuple layer thick Bi2Se3 thin film (mB = 0.32)(30j|. 
For a macroscopically large R, we found an approxi- 
mated solution for the energy spectrum of h+ as E = 
mjhvsgn(B) I R. As the time-reversal copy of h- has 
an approximated spectrum E = —mjhvsgn(B) / R. They 
form a series of paired helical edge states, in good agree- 
ment with the edge-state solutions in a 2D quantum spin 
Hall system 3l| if we take k — rrij/R for a large R. 
When shrinking R, the energy separation of these edge 
state AE = ±Hv /R increases with decreasing R, and 
the edge states with higher mj will be pushed out of 
the energy gap gradually. However, we observe that for 
mB > two pairs of states for rrij = ±1/2 always stay 
in the energy gap, and as R — > their energies approach 
to E = ±(v 2 /2|B|)V4mB - 1 for mB > 1/2 or ±mv 2 
for < mB < 1/2. The solutions demonstrate the for- 
mation of the in-gap bound states as illustrated in Fig. 
[TJ Therefore considering the symmetry between h + and 
h- we conclude that there always exist at least two de- 
generated pairs of bound states in the energy gap in 2D 
quantum spin Hall system in the presence of vacancy. 

The mechanism of the formation of the in-gap bound 
states is applicable to 3D topological insulators. 23, 32[ 
In 3D, the modified Dirac equation with a central po- 
tential becomes a classical problem, the hydrogen atom- 
like problem. For the Coulomb potential, it was exactly 
solved to give the fine structure of light spectra of hydro- 
gen atom. Similarly, the eigenstates of the 3D modified 
Dirac equation with a central potential can be labeled by 
three good quantum numbers. The first two are the total 
angular momentum J = fxp+-|S and its z-component 
J z , where the spin operator S Q = ao ® u a (a = x, y, z). 
The eigenvalues of J 2 and J z are j(j + l)h 2 and rrijh, 
respectively, with j g {h, §,•• •} and rrij G {—j, ■ ,j}- 
The third conserved quantity is the spin-orbit operator 
k = /3(f X p • £ + H). Note that k 2 = P + ^ 2 /4, then 
the eigenvalues of k is Hk — ±h(j + 1/2) = ±h, ±2h, 
Thus k here is similar to the ± index in 2D that sep- 
arate the hamiltonian into h±. These conserved quan- 
tities also help to reduce the problem into a set of ID 
radial equations. 33j In the presence of the vacancy or 
a cavity of radius R with the boundary conditions at 
^(R) = ^(oo) = 0, the radial part of the wave function 
can be solved in terms of the modified spherical Bessel 
function of the second kind k n (x). With the help of the 
recursion relation of k n (x) , the transcendental equations 
for the bound state energies can be found as 
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for k — j + i and — (j + ^), respectively. The correspond- 
ing wavefunction ^ . 3 K (r, 9, </>) are of the form 
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fc. T l 
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(9) 



(10) 



and Y7 



and 



■ , y 3, ip) is the spherical harmonics, 
possess opposite parities. 

Although the rotational symmetry simplifies the prob- 
lem, it is believed that the presence of the bound states 
is not sensitive to the shape of the vacancy, because of 
their topological origin. As an example we choose a set 
of parameters based on first principles calculations for 
Bi2Se3 by ignoring the anisotropy, where mv 2 — 0.28eV^, 
hv = 3.2eVA, and B = 33eVA 2 . In this case mB ~ 
1 > 1/2. Similar to the 2D case, we find that the sur- 
face states around the cavity exist for a large radius R 
as expected by the bulk-boundary correspondence for a 
Z2 invariant topological insulators (29j|. The states with 
larger orbital angular momentum are eventually expelled 
into the bulk band while the radius is shrinking. We 
plot several bound state energies of small orbital angular 
momenta as a function of the radius R in Fig. [3l For con- 
venience, the bound states are labeled by the quantum 
number k for the spin-orbit operator. Each k corresponds 
to (2j + l)-fold degenerate states of different rrij. Note 
that when the vacancy radius is only several angstroms, 
two degenerate pairs of bound-state energies can survive. 
Detailed analysis of the solution indicates that the spa- 
tial distribution of a bound state is comparable with that 
of the edge or surface states (for a large R in the present 
case), which are determined by the model parameters and 
slightly depends on R. From the evolution of the edge 
or surface states into the in-gap bound states, we think 
their formation have the same topological origin. Thus 
these in-gap bound states are also protected topologically 
as the edge or surface states. 

Now we come to address possible implication of these 
solutions to topological insulators. Due to the overlap- 
ping in energy, when the vacancies or defects are lo- 
cated close to the boundary, the induced in-gap bound 
states may sabotage the electronic transport through the 
boundary states. When the wave functions of the in-gap 
bound state and the edge or surface states overlaps in 
the space, the distortion of the wave functions of these 



4 



0.50 
0.25 
0.00 
-0.25 
-0.50 



— ^ 


3 


4 




- 


2 

K 1 






- 










-1 ^ 








-3 


-4 







I ' I ' 
1 2 


i 1 

3 






R (nm) 







negative effects of impurities and vacancies are precisely 
balanced. The topologically-protected bound states for 
sure are essentially different from those we know before, 
as they are subjected to some topological nature and con- 
fined to a mesoscopic scale. Their possible impact and 
applications for topological insulators in future deserve 
further studies to explore. 
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FIG. 3: Three-dimensional topologically protected bound 
states. Energies (E in units of the band gap A) of in- 
gap bound states covering a vacancy in a 3D topological 
insulator as functions of the vacancy radius R. k is the 
quantum number of the spin-orbit operator. Parameters: 
mv 2 = 0.28eV, hv = 3.2eVA, and Bti 2 = 33eVA 2 . mB = 0.90 
and A = 0.50eV. 



states due to the boundary conditions will cause the en- 
ergy change in these states. As a result, we may also 
regard that there exists a transition amplitude between 
these states. For example, in the 2D quantum spin Hall 
system, there exists a pair of helical edge states with a 
linear dispersion, Ek.± — ±hvk. Consider a vacancy or 
defect appears near the boundary. The effective model 
for a pair of helical edge states in the presence of the 
defect states has the form, 



H 



hvkacl^Cka- 



k.a.n 



(T^cldn+h.c.) 



where c ka and c^a are the creation and annihilation oper- 
ators of edge states and d^ n and d n are for the nth in-gap 
bound states with the energy e„. The resulting disper- 
sions for the edge states are no longer linear in the mo- 
mentum, and open energy gaps AE = 2 \T£ a \ at the res- 
onant point e n = hvka. TJ} a is a function of the relative 
position between the defect or vacancy and the boundary. 
In the HgTe/CdTe quantum wells, a typical size of the 
edge states and bound states is about 50 - 100 nm. The 
energy gap opening here is not caused by breaking the 
time reversal symmetry. The mechanism is similar to the 
finite size effect in the quantum spin Hall system, where 
the overlap of the edge states living on opposite edges 



causes an energy gap[3lj. In that case the energy split- 
ting is about 0.5 meV for a strip system with a width of 
200 nm. The effect is large enough to be measured exper- 
imentally. This may help to understand why the non-zero 
conductance is narrowed to a small region of gate voltage 
in the HgTe/CdTe quantum wells [22]. The in-gap bound 
states may also be one of the possible mechanisms for the 
low mobility in 3D topological insulators. 

However, blessings usually come in disguise. The whole 
semiconductor business depends on how the positive and 
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